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ABSTRACT 


Atmospheric optical turbulence induces phase fluctuations in a propagating elec- 
tromagnetic wave. The resulting degradation in coherence limits the capability of any la- 
ser, target acquisition, or surveillance system. Past data collection methods for the 
parameterization of atmospheric turbulence profiles, in support of critical Theater Ballistic 
Missile Defense (TBM) systems, from ground level to 30 km, have depended on mete- 
orological balloon-thermosonde systems, probes carried on the U.S. Aur Force Argus air- 
craft, as well as radar and optical measurements.. The balloon and aircraft systems 
measure temperature fluctuations and compute the temperature structure function, Cr’ and 
the related index of refraction structure parameter, C,”. It has recently become critical to 
explain why turbulence profiles from daytime thermosonde data consistently show a two 
order of magnitude increase over that taken during the night, primarily between 12-20 
km. 

This thesis analyzed the TSI 3.8 um platinum coated tungsten thermosonde probe 
used by the USAF Research Laboratory (AFRL) to quantify the magnitude of the solar 
heating and to investigate other heat transfer mechanisms in the probe. A model of the thin 
wire probe was developed to identify each of the contributions to the temperature error 
and its significance. Experimental measurements where collected to verify most aspects of 
the final model. 

We found that the sun induces a temperature mse in the TSI 3.8 um fine wire 
probe, during the day, that can vary from near zero to 0.175 K. It is strongly dependent on 
probe orientation with respect to the sun and on variations in the air flow over the probe. 
This then causes an apparent increase by two orders of magnitude in the daytime meas- 


urements of the optical turbulence parameters Cr’ and C,/. 
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a: INTRODUCTION 


The Airborne Laser Program (ABL) must contend with several problems associ- 
ated with the propagation of electromagnetic radiation through the atmosphere. Specifi- 
cally, random fluctuations in the atmospheric index of refraction cause phase front 
distortions of the propagating plane wave [Ref. 1:p. 6-6]. Other effects such as thermal 
blooming, absorption and scattering by atmospheric particles, and stimulated Raman scat- 
tering are also important but go beyond the scope of this thesis [Ref. 2]. The concept of 
mounting a weapons class laser on an aircraft has been around for over two decades and 
was proven possible by the Airborne Laser Laboratory in the 1970’s. The technologies to 
make it practical, however, have only recently been demonstrated by Lockheed Martin and 
scale laser testing done by the Air Force Research Laboratory in 1995 as part of the 
ABLE ACE experiment [Ref. 3:p. 6]. 

The basic requirements of the ABL system are to detect, track, and destroy theater 
ballistic missiles (TBMs) while still in their boost phase and prior to the missile’s leaving 
the hostile country’s own airspace. The ABL system will use a passive infrared (IR) 
search and track system along with an active Multi-Beam Illuminator (MBIL) laser for 
accurate positioning and better dwell characteristics during high power Chemical Oxygen 
Iodine Laser (COIL) missile destruction [Ref. 4:p. 1]. Adaptive optics will be needed to 
achieve effective ranges beyond ~100km. 

The ABL program utilizes atmospheric turbulence data that have been collected 
by the Air Force Research Laboratory (AFRL), Naval Postgraduate School (NPS), and 
Argus aircraft for this and other Department of Defense atmospheric laser programs over 
the past decade. Characterization of atmospheric turbulence is essential in the design and 
use of atmospheric laser systems. A critical factor in the design is whether a perceived 
enhancement in the measured temperature structure parameter, Cy’, during the daytime as 
compared to night is real or is a consequence of solar heating of the measurement appa- 
ratus itself [Ref. 5:p. 30]. Experiments done by AFRL, as well as subsequent ones done 


by NPS, have shown a significant enhancement in daytime Cr” values, most likely caused 


solar heating of the thermosonde probes used to do the measurement, coupled with the 
complications of package rotation, pendulum motion, and probe vibration [Ref. 5:p. 55]. 

This thesis investigated the systematic effects caused by solar heating of the bal- 
loon thermosonde probes used by AFRL and NPS. Extensive probe response modeling 
and experimental measurements were done to identify and to quantify the individual con- 
tributions to temperature variation and to determine if the signal generated would be in- 
terpreted as a false turbulence signal within the electronic bandpass of the instrumentation. 
Package rotation, pendulum motion, and other dynamics where apparent in the exper- 
mental data. 

The ABL program has entered the Program Definition and Risk Reduction phase 
(PDRR) which calls for one complete weapon system, including a laser and its controlling 
system, to be integrated onboard a 747-400 platform and to be successfully demonstrated 
by 2002. Success of the system will be evaluated by attempting to shoot down a SCUD 
type missile [Ref. 6:p. 1-2]. Essential msk assessment and management issues underlying 
the feasibility of proceeding past the project’s scheduled acquisition program Authority to 
Proceed Milestone (ATP) in mid-1998 relies partly on the resolution of the atmospheric 
turbulence problem. The underlying motivation of this investigation was to identify any 
systematic error sources that might affect the microthermal platinum coated tungsten 
probes that provide critical turbulence measurements needed for effective deployment of 


theater ballistic missile defense systems. 


I. BACKGROUND 


A. OPTICAL TURBULENCE PARAMETERS 


Atmospheric turbulence degrades the wave front of a collimated laser beam. Ran- 
dom fluctuations in the atmospheric index of refraction, n, caused by temperature and 
humidity inhomogeneities, induce centroid wander, beam broadening, and interference 
that all reduce the beam’s spatial coherence (Figure 2-1) [Ref. 7:p. 1]. 

Tatarski and others have developed useful parameters for the statistical quantifi- 
cation of optical turbulence in the atmosphere by defining so-called structure functions. 
By describing the atmosphere as a three dimensional random field Tatarski considered the 


correlation function of some relevant quantity of interest, f, as [Ref. 7:p. 3] 


Bett) =(FG)-(Fq)) 7° G)- (4° @))), (1) 
where the angled brackets denote the mean value, the superscript star 1s the complex con- 


jugate, and t denotes time. Because the atmosphere 1s not statistically homogenous or 
stationary, the mean values ( Ff (7, ) are not well defined. To avoid this problem Kol- 
mogorov and Tatarski introduced the structure function 

Dy (t,.t2)= ([F,)-F@)] ). (2) 
which characterizes the intensity of the particular random fluctuations of f(t) [Ref. 7:p. 9- 


10]. For our use then, a structure function tensor can be described by the mean square 


difference between two measurements, such as temperature, by 


D,(%%) = (TG) - THY). (3) 


Initial Plane Wave Turbulence Pockets Resultant Distortion of Plane Wave 





Figure 2-1: Effect of atmospheric turbulence on a coherent plane wave. 


For an isotropic and homogeneous quantity the structure function must necessarily be in- 
dependent of the onentation angle and absolute position, and can depend only on r, the 


distance between the two observation points: 

r=|R-Al, (4) 
Although the atmosphere is generally not isotropic or homogenous, if r is sufficiently 
small, then some local volume may be sufficiently isotropic and homogenous to permit 
credible measurements [Ref. 7:p. 31]. 

It is found that the statistical behavior of the atmosphere is well behaved over a 
limited range of r, 1, <r < Z,. This is called the inertial subrange. The inertial subrange 
is bounded on the lower end by /,, which is the small scale turbulence region associated 
with viscous eddy dissipation. On the upper end of the inertial subrange, Lo defines the 
outer scale length, which is the largest eddy size that can be considered isotropic through- 
out its volume [Ref. 8]. 

Kolmogorov and Obukhov proposed that the structure functions in the atmosphere 
have a r~> dependence [Ref. 7:p. 31]. Assuming isotropic and homogenous conditions, 


4 


Tatarski proposed the temperature structure function, Cy’ by normalizing Eq. (3) by r° 
[Ref. 7:p. 46] 


co - Dil) _ (7G)- TEP) . 


i is 





wi 


Since 7-1 1s proportional to density, which is a function of the pressure, P, in 
millibars and temperature, T, in K, the inertial range refractive index structure parameter, 
C,”, is related to the temperature structure parameter by 


P(z) 
T’(z) 





Gag) = 799210 | (2), (6) 


where z is the altitude in meters [Ref. 7:p. 79]. The temperature structure parameter, Cy” , 
is measured by fine wire temperature sensors in a balloon or aircraft system as a function 
of altitude. C,” is then computed from Cy” using Eq. (6). 

The effects of turbulence on a laser beam propagating through the atmosphere ap- 
pear as distortion in the final beam pattern at the target. A common measure of beam 
quality is the Strehl ratio, which is the relative intensity of the beam’s point spread func- 
tion measured at the mean centroid location relative to that of an ideal, diffraction limited 
optical system [Ref. 9]. The Strehl ratio depends on the spatial coherence of the laser ra- 
diation after passing through the atmosphere, which in turn depends on the atmospheric 
turbulence and its effects on the modulation transfer function [Ref. 18]. The atmospheric 
modulation transfer function (MTF), characterized by the transverse coherent path length 
Io, is a measure of electric field transverse autocorrelation length, [Ref. i:p. 6-29]. Fried 


showed that, for plane waves, 
z -3 
r, = 21146k? | C2(z)az]? , (7) 
0 


where k is the optical wave number (27/() and L is the optical path length [Ref. 10:p. 


397]. Values of ro range from centimeters under very turbulent conditions to meters 





Infinite coherence length Turbulence Pockets Resultant reduction of ro 


Figure 2-2: Physical representation of the coherence length, ro, and its reduction due to 
atmospheric turbulence. 

when turbulence is very low. Since ro describes the spatial coherence of the beam, the 
larger its value the better (Figure 2-2). The diffraction limit of the system, however, de- 
pends on the geometry of the optical system as well as ro. Without turbulence the diffrac- 


tion-limited resolution in radians for a circular aperture is [Ref. 11:p. 393] 


7) 


ae ere (8) 


where A is the wavelength and D is the aperture diameter. When D>ro, the effective tur- 


bulence-limited angular resolution just becomes 


A@,in = —- (9) 


Lo 
The isoplanatic angle, 80, is another useful statistical parameter for the characteri- 


zation of atmospheric optical turbulence. It is related to mr, and represents a measure of 


the angular field of view of an optical system. It characterizes a cone within the atmos- 


6 


phere with similar turbulence characteristics. Fried [Ref. 12] showed that 6) depends on 


the integration of C,” along an optical path and is given by 
L -3/5 
0, = 29050 | Cc Ode (10) 
0 


with L, the optical path length in meters. 

As has been previously stated, current balloon and aircraft probe thermosonde 
systems can measure Cy’ and then compute C,’ as a function of altitude and time. Using 
Eqs. (8) and (9) this data provides the basis for computing ro and Qp for a particular appli- 


cation in the design and development of an adaptive optics system. 
B. THERMOSONDE 


The AFRL thermosonde design used during the last decade consists of a VIZ 8000 
electronics package with two fine wire TSI 3.8 tum platinum coated tungsten probes ex- 
posed to the ambient atmosphere and separated by 1 m (Figure 2-3). This system meas- 
ures the wires’ resistance fluctuations caused by changes in temperature. The balloon 
carrying the package is capable of ascending to an altitude of approximately 30 km. 

Changes in the probe resistance, Rw, give the temperature, Tw, through the rela- 
tion 

R, = R, +[1+ a(Z, -%)] (11) 
[Ref. 13:p. 11], where Ro and To, are values at a specific reference resistance and tem- 
perature respectively and ao is a constant for the particular wire composition and geome- 
try, platinum coated tungsten in this instance. In practice, the two probes form the 
opposing elements of an AC Wheatstone bridge. The difference in resistance from the 
Wheatstone bridge is amplified and converted to a temperature change by [Ref. 14:p. 15] 
_(R,-R) 


AT oR, 


(12) 


The probe data channels are sampled every 10 s along with ambient pressure, mean tem- 
perature, and humidity. The data are radioed via the modified VIZ 8000 radiosonde to the 
ground station that continuously displays and stores the data on a PC. Optimum ther- 
mosonde performance requires a compromise in the probe current. A value of 0.1 ma for 
a 20 Q probe provides sufficient signal without inducing temperature differences above 
0.002 K between the probe and the atmosphere. The ohmic heating temperature differ- 


ence, AZ, is 


q., 
AT, =~, (13) 





where qyw is the electrical power input, A, is the surface area of the wire, and h is the coef- 
ficient of heat transfer (Ref. 14:p. 17]. 

The response of the probe to a change in temperature must be faster than the ac- 
tual fluctuation to be measured if a meaningful measurement is to result. Failure of the 
probe response to follow the actual atmosphere will cause an underestimate of the struc- 
ture function. The response time of our probe is given as [Ref. 14:p. 20] 


C 
4 pDe 


= ay an as 





with: 


p = density of probe wire (kg/m’), 

c= heat capacity of wire (J/kg K, 

h= convective heat transfer coefficient (W/m’K), 
D= wire diameter (m). 
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Figure 2-3: Construction layout of platinum coated tungsten AFRL probe and support. 


G. CURRENT APPLICATIONS 


The apparent increase in turbulence during the day compared to night has signifi- 
cant implications on the ABL program (Appendix A). Current performance specifications 
for the ABL’s COIL laser require the destruction of a ballistic missile that is still in it’s 
boost phase out to a range of several hundred kilometers [Ref. 15:p. 1]. The power re- 
quirements on target for these ranges have only recently become possible with the use of 
adaptive optical systems that rely on active compensation of wave front distortion of the 
propagating laser beam. Still, only under favorable atmospheric conditions can ranges 
such as these be attained. If existing daytime balloon measurements of Cy’ are valid they 
imply that the useful range of the laser system is severely limited in the daytime. Other new 
programs which rely on laser propagation through the atmosphere to great distances 


would also have to be substantially re-thought. 


~ 


Dy: PREVIOUS RESULTS 


Numerous experiments have been performed by AFRL and NPS to characterize 
the microthermal probes and to determine Cy’ and C,” during day and night conditions. 
Brown has flown many balloon launches since July of 1981 and performed extensive 
analysis to show that the probe supports do not contribute to observed diurnal differences 
[Ref. 5, 14, and 16]. Walters, Weitekamp, Roper, and others through multiple balloon 
launch experiments (notably from 31 May to 3 June 1990 at AMOS observatory on Mt. 
Haleakala, Hawan and 4 -9 December 1990 at Holloman Air Force Base, New Mexico) 
showed that solar contamination of data was a real effect. Figure 2-4 is a representative 
example of night time data collected over Maui on 2 June 1990 showing that at least some 
of the raw C,” data reached the lower noise limit (solid line) of the sonde electronics. Fig- 
ure 2-5, shows a two order of magnitude shift in C,” data for a next day launch from the 
same location. None of the data approached the noise limit of the electronics. One would 
expect that under normal conditions, there should be some parcel of air within the atmos- 
phere that is not turbulent. Yet the measurements show that that assumption is valid at 
night but not during the day. Hence some doubt has been cast on the general validity of 
the daytime measurements. 

Additional work by Walters, Weitekamp, and Roper have also shown that balloon 
wake effects are not a source of contamination [Ref. 17 and 18]. Additional attempts at 
methods of sonde stabilization have also been investigated by them to eliminate solar ef- 


fects, but to date no satisfactory method has been found. 
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Figure 2-4: Raw night time AFRL probe C,. data collected on Maui, 2 June 1990 show- 


ing the lower electronic noise limit (solid line) and corresponding data points near this. 
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Figure 2-5: Raw day time AFRL probe C,’ data collected on Maui, 2 June 1990 showing 
the lower electronic noise limit (solid line) and corresponding data points with solar in- 
duced shift of two orders of magnitude. 





Hi. MODELING THERMOSONDE PROBE HEAT TRANSFER EFFECTS 


A. INTRODUCTION 


In this Chapter, the heat balance equation for the TSI probe wire is described in 
detail. Special attention is paid to convective cooling of the wire, the effects of changing 
solar illumination areas, and effective free stream velocity fluctuations. MATLAB pro- 
gramming was used to model the probe under these changing conditions as a function of 
altitude. Once the model was validated by experimental measurement as described in 
Chapter IV, simulations of different probes’ responses to material and geometry changes 


was investigated as shown in Chapter V. 


B. HEAT BALANCE EQUATION 


The equilibrium temperature of a fine wire probe is given by the heat balance 


equation 

OQ sonv, take ter > Cree me oye — 0, (15) 
where: 
Qconv = Convective transfer of power from the probe to the atmosphere (W), 


Qsoiar = Solar power absorbed directly by the probe (W), 

O23 reflected solar power from the Earth absorbed by the probe (W), 

Qpianck = Planck radiation power transfer between the probe and the environment (W), 
Qohmic = Ohmic heating caused by current in the probe wire (W). 


I 


In terms of power (W) the equilibrium condition is 


—RA, AT + 0A, P, +0A,P.y +(€A-P., - €0A,T* + €A4sP,)+?R=0. (16) 
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The individual terms in Eq. (16) will be defined in the section in which the heat 
transfer mechanism is introduced. The thin wire probe areas A; -As depend on the heat 
transfer mechanism, package pendulum motion and rotation, probe wire vibration, and 
solar zenith angle and are also defined in their appropriate sections. 

The time for wire to reach thermodynamic equilibnum is short compared to the 
heat transfer time via convection, so equilibrium was assumed [Ref. 17:p. 32]. Thermal 
conduction through the probe wire to the supports was also ignored as shown in [Ref. 
S:p. 45] for the simulation, but was experimentally measured as described in Chapter IV. 

Dunng flight, the effective free stream air velocity over the probe wire is also vari- 
able. Variations in the wind velocity alter the heat transfer from a probe that has a different 
temperature from the air. These fluctuations were modeled both as a steady state charac- 
teristic air velocity (V= 5 m/s) and as a random fluctuation distribution about a character- 


istic free stream velocity (Figure 3-1). 
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Figure 3-1: Simulation of air velocity fluctuation over probe 
wire for a fluctuation of 2.5 m/s about a velocity of 5 m/s. 
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ce. CONVECTIVE COOLING 


1. Convection with Steady State Free Stream Velocity 


The convective transfer of power from the probe wire to the atmosphere is 


Qeonv=hA AT, ( 7) 


where: 


h = convective heat transfer coefficient (W/m K), 
A, = aDL = surface area of the probe subject to convection (m’), 


AT = ne i i (K). 


Heat transfer from a heated probe wire to the atmosphere consists of free convection from 
buoyancy effects at low air velocity, forced convection at higher air velocities (Reynolds 
numbers, Re >0.5), and mixed free and forced convection between both these regions, 
about Re=0.5 [Ref. 20:p. 323]. Figure 3-2, developed by Collis and Williams, illustrates 
the complicated behavior of this low Re region. 

The Reynolds number is the ratio of inertial to viscous forces on a parcel of fluid. 


In this application, 
Re= sae (18) 
Vv 


where D is the probe diameter, v, is the kinematic viscosity (which is the ratio of dy- 
namic viscosity, 1, and density, p) and V is the characteristic velocity [Ref. 21:p.54]. The 
NPS-AFRL probe, as modeled, would experience Reynolds number as low as 0.00024 but 
no more than about 4. The Nusselt number 


hD 
Nu = a (19) 


t 


is another useful dimensionless number describing the heat transfer coefficient, where k, is 


the thermal conductivity of the fluid, in this case air. 
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Figure 3-2: Interaction of free and forced convection from 
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At low Reynolds numbers, where buoyancy becomes important, the dimensionless 


Grassof number, given by 


aie PEP Ten pe | o 
L 

relates the ratio of buoyant to viscous forces near the wire where g is the local gravita- 
tional acceleration and B2 is the coefficient of expansion of the fluid. The Grassof number 
is important when Re<Gr’” and buoyant effects become a significant contributor to the 
probe’s ability to dissipate heat [Ref. 22:p. 131]. Grassof numbers modeled for our probe 
varied from near zero to 3.25x10°. The effect of the Grassof number on the buoyancy ef- 
fect for low Reynolds number 1s also shown in Figure 3-2. At altitudes above 30 km mo- 


lecular diffusion becomes dominant over convective effects for the heat transfer 


mechanism of the AFRL probe wire [Ref. 14:p. 19]. 


2. Determination of h 


The convective heat transfer coefficient h is the heat flux leaving the wire surface 
per unit temperature difference between the wire and the free stream. The effect of h on 
the probe is significant and is a complex function of geometry, Reynolds, Grassof, and 
Prandtl numbers as well as the overheat ratio, (Tprote-T)/T [Ref. 23: p.352]. The dimen- 
sionless Prandtl number, Pr, is the ratio of the momentum diffusivity to the heat diffusivity 


and is given by 





i = , (21) 


The conductivity k, is defined in Table 3-1, and c, is the specific heat of the atmosphere. 
The convective heat transfer coefficient can be found using a variety of methods. 
Most relate the Nusselt , Reynolds, and Prandl numbers and are valid over certain Rey- 
nolds and Nusselt number regimes, probe geometries and fluid/gas compositions. No sin- 
gle credible model for our probe’s entire operating environment was found. Four different 


models were compared, and the relative merits of each were contrasted to our probe’s 
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actual performance. Consideration was also given to the physical processes that were in- 


cluded in the model. The four methods considered are summarized below. 


a. h, Using Kothandaraman and Subramanyan’s Method 


Kothandaraman and Subramanyan in [Ref. 24: p. 91] proposed the follow- 
ing relation for the Nusselt number 
hD 


Nu, = ae CRe’”, (22) 


t 


where k, is the conductivity of the fluid (air) determined empirically in Table 3-1: 


Alp 505i =) 0S" Alt + 206.5 


k. (W/mk) = 1983.47 + 11636 





Table 3-1: Conductivity k, of atmosphere variation with altitude in kilometers, adapted 
from the Geophysics Handbook [Ref. 25:p. 14-6]. 


C and n are constants as given in [Ref. 24:p. 91] for the Reynolds number range 
0.4<Re<4.0 and have values of C=0.891 and n=0.330. Values of Re less than 0.4 require 
further extrapolation, as done in [Ref. 17]. The altitude, Alt, here and throughout, is 
given in kilometers. A characteristic velocity V of 5 m/s, corresponding to the approxi- 
mate ascent rate of the balloon, and characteristic length D of 4.7 um (measured fine wire 
probe diameter; refer to Chapter V) were used. A more elaborate treatment for modeling 
the effective velocity over the probe wire was done in Ref. 5:p. 7-9 and Ref. 19:p. 22, but 
for our methods this characteristic velocity will be sufficient. The density assumed the va- 
lidity of the ideal gas law 

PRT 


= ——— 23 
va (23) 
and modeled the pressure variation with altitude by 
P= pe*. (24) 


18 


R, above, is the universal gas constant (8.31 J/Mole-K) and M is the mean molecular 
weight of the atmosphere (28.9644 kg/kMol) [Ref. 25:p. 14-5]. Atmospheric tempera- 
ture variation with altitude, T, came from the USAF Geophysics Handbook and is shown 
in Table 3-2. The scale height, z, and scale pressure, po, are 8.1 km and 1013 mBar 
(101.3 kPa) respectively [Ref. 25:p. 14-5]. 
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Table 3-2: Atmospheric mean temperature variation with altitude adapted from the Geo- 
physics Handbook [Ref. 25:p. 14-3]. 


Dynamic viscosity was determined using the relation 
El 
AL? 
TG 
[Ref. 25:p. 14-5]. Here B; is a constant (1.45x10° kg/s-m-K”’) and S is Sutherland’s 
constant (110.4 K). 





(25) 


With all quantities known and solving for h we find 
CRe*k 
D 


a 





(26) 
For this probe the diameter, D, was 4.7 um and the length, L, was 4.15 mm. 


b. h, Using Kramer’s Method 
Kramer’s law for circular cylinders of infinite length is given by Hinze in 
[Ref. 26:p. 201] as 
Nu = 0.42 Pr°” +057 Pr°* Re®”’. (27) 
Kramer’s law is valid for Reynolds numbers in the range 0.1<Re<10,000 [Ref. 19:p. 63]. 
So now using Eq. (23), h can be found. 


= 
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c: h, Using Zijnen’s Method 


The heat transfer proposed by Zijnen does not have any stated valid Rey- 
nolds number range for circular cylinders and is given as [Ref. 19:p. 64] 


Nu = 0.38Pr°? + (056Re° + 0.001 Re) Pr? (28) 


Once again Eq. (23) is used to solve for h. 


d. h, Using Collis and William’s Method 


Collis and Williams developed the following relation for the Nusselt num- 
ber [Ref. 27:p. 357-384] 


(29) 


T A a 0.17 


Nu, =(A ~BRe}) 
For air and in the Reynolds number range 0.02<Re<44 the constants are as follows: 
n=0.45, A=0.24, and B=0.56 [Ref. 22:p. 130]. The subscnpt f denotes measurements 
were taken at the film temperature which is the average temperature between the wire and 
the atmosphere. Eq. (23) then gives the relation for h in this method. 

Baille further modified this relation to account for the mixed convection effects in 
the region around Re<0.15, as shown in Figure 3-2, by introducing an effective Reynolds 
number, Reg, given by 

Rew — Re occr =)", (30) 
Baille also investigated the error associated with assuming an infinite cylindrical probe in 
the above approximations and developed Figure 3-3, which shows the effect of aspect ra- 
tio (21/D) on the Nusselt number [Ref. 22:p. 132]. The AFRL probe used has an aspect 


ratio of over 2000 and therefore behaves as an infinite wire to high degree. 
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Figure 3-3: Effect of wire aspect ratio on Nusselt 
number. From [Ref. 22:p. 132]. 


e. Comparison of h Values 


Although no method covered the extremely low Reynolds numbers en- 


countered in this work, the Collis and Williams method best represented the physical ef- 
fects encountered and the behavior of our probe that is as most represented by the low 
Reynolds number regions shown 1n Figure 3-2. All four methods were in good agreement 
at Reynolds numbers greater than 1. Final values for the coefficient of heat transfer for air 


velocities of 0 and 5 m/s as derived from the Collis and Williams technique are shown in 


a 


The kink in the curves in Figures 3-4 to 3-9 around 12 km comes from the rate of 
change of temperature versus altitude of this part of the stratosphere. 

Figure 3-5 summarizes the coefficient of heat transfer values for the four methods 
as a function of altitude for an air velocity of 5 m/s using a MATLAB program (Appendix 
B). This program was also used to produce the plots of the, Re, Nu, and Gr numbers, at- 
mospheric temperature T and kinematic viscosity v versus altitude, as shown as Figures 3- 
6 to 3-11 for the steady state free stream velocity of 5 m/s and as a function of altitude. 

Although not a factor with balloon mounted probes, those mounted on aircraft, 
such as Argus, also need to consider the heat loss effects at different Mach numbers. 


Spangeuberg [Ref. 22:p. 133] presents the reported heat loss at various Mach numbers 


and is reproduced here as Figure 3-12. 
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Figure 3-4: Convective heat transfer coefficient of 4.7 um probe 
using Collis and Williams method versus altitude with constant 
free stream velocity of 0 and 5 m/s. 
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Figure 3-5: Comparison of h, for the 4.7 um probe, using 
the four methods investigated at an air velocity of 5 m/s. 
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Figure 3-6: Reynolds number versus altitude, for the 4.7 um probe, 
for a constant free stream velocity of 5 m/s. 
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Figure 3-7: Reynolds number versus altitude, for the 4.7 um 
probe, for an air velocity of 0.01 m/s. 
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Figure 3-8: Nusselt number versus altitude, for the 4.7 um probe, 
for a constant free stream velocity of 5 m/s. 
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Figure 3-9: Grassof number, showing ratio of buoyancy to vis- 
cous forces, versus altitude for constant free stream velocity of 
5 m/s and corresponding to a differential probe temperature of 


0.035 K for the 4.7 um probe. 
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Figure 3-10: Atmospheric temperature versus altitude. 
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Figure 3-11: Kinematic viscosity of the atmosphere versus altitude. 
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Figure 3-12: Heat loss at increasing Mach number, from Davies and 
Fischer in [Ref. 22:p. 133]. 
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2. Convective Cooling Power 
Typical values for the convective cooling of the wire at sea level are on the order 
of 10.63x10° W, but they are sensitive to the free stream velocity. Shown as a function of 


altitude in Figure 3-13 are values of Qconv for a free stream velocity of 5m/s. 
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Figure 3-13: Probe convection, Qconv versus altitude (V = 5 m/s) 


4, Convection with Normal Free Stream Velocity Fluctuation 


Real world conditions encountered by the probe involve dynamic fluctuations of 
the free stream velocity component over the probe wire. These fluctuations arise from 
package pendulum motion, rotation, possible balloon wake turbulence, actual atmospheric 
turbulence, and to a lesser extent probe vibration. 

Velocity fluctuations ultimately affect the coefficient of heat transfer and Reynolds 
number which both depend on air velocity over the probe wire. The effect of a random 
normal mean free stream velocity fluctuation of 2 m/s about 5 m/s 1s shown in Figures 3- 

14 and 3-15 for both these quantities. The convection, Qconv, term will vary similarly with 
altitude and is not plotted here. 
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Figure 3-14: Modeled convective heat transfer coefficient with 
velocity variance of 2 m/s about 5 m/s for the 4.7 um probe. 
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Figure 3-15: Modeled Reynolds number with normal mean ve- 
locity variance of 2 m/s about 5 m/s for the 4.7 um probe. 
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D. DIRECT SOLAR ABSORPTION OF ENERGY 


1. Constant Probe Cross Sectional Area Presented to Solar Radiation 


The power absorbed directly from the sun by the probe wire is given by 

Qsotar=HA2P,, (31) 

in the heat balance equation with: 

a= solar absorptivity of the probe, 

A= cross sectional area of probe seen by the sun (m’), 

P,= solar power flux (W/m’). 

The solar absorptivity of the platinum coated tungsten wire probe is 0.435 from tabulated 

values of the normal incident reflectivity (R=0.565 @ 3.0eV) for platinum wire, which is 

related to the absorptivity by a=1-R [Ref. 28:p. 10-282]. The maximum illuminated area 

of the probe as seen by the sun to good approximation is the cross sectional area DL. The 

projected area, A>, is a function of the complex probe pendulum motion and rotation, 

probe vibration, and solar zenith angle. The change in area subjected to solar radiation, as 

a function of time, was modeled as one half the magnitude of the sum of two sinusoidal 

variations of the total visible surface area. Periods for the two sine functions correspond 

to the observed frequency of pendulum motion and rotation. Theta dependence (sun ze- 

nith angle) was not modeled explicitly in this thesis but was done previously [Ref. 5:p. 10]. 


The solar flux, P, varies with altitude (meters) as given in [Ref. 25] 
P =(Alt)’"" +890, (32) 


to our maximum altitude of interest, 30 km. 
Typical values of the power absorbed directly from the sun by the probe at sea 


level are of the order 7.28x10° W and are shown as a function of altitude in Figure 3-16. 


ZO 


30 


29 


Altitude (km) 
— N 
On Oo 


=-s 
Oo 


7 7-0 8 8.5 9 9.5 10 10.5 11 
Direct Solar Absorption, Qsolar (W) 


Figure 3-16: Direct solar absorption, Qyolar, Of power with constant 
exposed cross sectional area versus altitude for the 4.7 um probe. 


2. Variation of Probe Cross Sectional Area Presented to Solar Radiation 


The probe cross sectional area exposed to solar illumination is a function of the 
roughly sinusoidal geometry variations caused by pendulum motion and rotation with re- 
spect to the sun. Secondary influences of probe wire vibration are discussed below. Us- 
ing MATLAB programs developed by Walters (Appendix C), digital filters and Founer 
spectral techniques were used to analyze single probe thermistor data taken at Anderson 
Peak in September 1992. This launch was during the day and showed the characteristic 
C,’ enhancement obtained from higher order structure functions (Figure 3-17). Figure 3- 
18 shows the resultant power spectrum plot with two distinct periodicities of about 5 and 
20 seconds. Visual observation of these balloon launches and simple calculation of the 
pendulum motion confirm that these correspond to the periods of the rotational and pen- 


dulum motion of the sonde package during its ascent. Hence these were used in our 


- 


30 


model. The influence of the sinusoidal variation of the solar cross sectional area A> on 


solar absorption Q.olar is Shown in Figure 3-19. 
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Figure 3-17: Anderson Peak C,” during the day obtained 
from a single flake thermistor. 
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Figure 3-18: Power spectral plot of Anderson Peak data 
- from Figure 3-17. 
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Figure 3-19: Effect of Sinusoidal variation of probe solar illumina- 
tion area, A> on solar heat flux with air velocity of 5 m/s for the 4.7 
iim probe. 


S: Probe Vibration 
Less well understood 1s the effect on solar heating caused by the probe wire vi- 
brating in and out of solar view. Experimental work done to measure this phenomenon is 


discussed further in Chapter IV. 
E. REFLECTED SOLAR ABSORPTION 


The reflected solar power from the Earth onto the probe wire 1s given by 
Qrer=XAUP sy, (33) 
where: 


a= solar absorptivity of the probe, 

A, = effective area of the probe that absorbes reflected sunlight (m’), 
P.= solar energy flux (W/m’), 

Y= percent of solar energy reflected by the Earth (terrestrial albedo) 


Bz 


The solar absorptivity and solar flux terms are as previously determined. Solar flux is 
taken to be a fraction of that at the surface of the Earth (890 W/m’). The maximum cross 
sectional area of the probe seen by the earth is DL. The percent of solar power reflected 
by the Earth (albedo) was estimated to be 39%. 

Typical values of the reflected power of the sun from the Earth on the probe are of 


the order 2.84x10° W at sea level and are approximately constant with altitude. 


iF PLANCK RADIATION 


i Planck Radiation from the Atmosphere 


Planck radiation consists of three terms, that being absorption from the surround- 
ing atmosphere, that being absorbed from the Earth, and that being radiated away by the 


probe. It is represented by 


Qolanck=QatE-Qp. (34) 
Energy transfer in the form of Planck radiation from the atmosphere to the probe is 
qa=EAcP atm, (3 5) 


where: 


€= emissivity (0.065), 
Ac= area of the probe presented to the atmosphere (m”), 
Piim= thermal radiation emitted from the atmosphere (W/m’). 


The emissivity term is taken as a constant and given the value of that for 300 K radiation. 
The area Ag is the total probe surface area. The radiation term Pr, was estimated by inte- 
grating over the plot of radiance versus wavelength for sky radiation as seen from the 
Earth, given in [Ref. 29:p. 3-75]. For a first order approximation the atmospheric flux on 
a sunlit summer day at 30°-40° latitude with an estimated probe field of view of 47 stera- 


dians was used. An approximate value for Pam of 250 W/m’ resulted. 


Si 


Typical values of radiation from the surrounding atmosphere on the probe are of 
the order 9.6x10°’ Wat sea level and are constant as a function of altitude. Pam is more 
correctly a function of the air temperature and therefore would vary with altitude, but 


these effects are small and were neglected in this work.. 


o Planck Radiation from Earth 


Planck radiation from the Earth is given by 
qe=eAsP., (36) 
with: 


€= emissivity for 300K radiation, 
As= area of probe presented to earth (m’), 
P= radiation from Earth (W/m’). 


Emissivity is as before and As is equal to the cross sectional area DL. The Earth-radiated 
flux was estimated using the same integration approximation technique as used for Patm 
above but now integrating over the upper portion of the radiance versus wavelength plot 
of the Earth as seen from space as given in [Ref. 29:p. 3-54]. This integration yielded a 
first order approximation of the solar flux on a sunlit summer day at 30°-40° latitude with 
an estimated probe field of view of 27 steradians. An approximate value for P. of 615 
W/m’ was used. 

Typical values of radiation on the probe are of the order 7.5x10°’ W at sea level 


and are constant as a function of altitude. 


3. Planck Radiation of the Probe Wire 


Planck radiation from the surface of the wire is 
qp=ecA3T", (37) 
where: 


E= emissivity for 300K radiation, 

o= Stefan-Boltzman constant (5.67x10° W/m°K’), 
A3= area of the wire probe radiating (m’), 

T= atmospheric temperature (K). 
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The emissivity for platinum wire was 0.065 [Ref. 28:p. 10-282]. The area term, As, is the 
total surface area of the probe. The atmospheric temperature is used instead of the probe 
temperature Tprote to simplify the calculation and because AT<<Tp,obe aS Shown in [Ref. 
S:p. 4]. 

Typical values of the Planck radiation emitted by the probe are of the order 


1.50x10° W at sea level and are shown as a function of altitude in Figure 3-20. 


30 
25 


20 


Altitude (km) 
or 


10 
5 
0 
0.4 0.6 0.8 1 1EZ 1.4 1.6 
Planck Radiation by the Probe Wire (W) x 10° 


Figure 3-20: Planck radiation from the 4.7 um probe as a function 
of altitude with constant air velocity of 5 m/s. 


G. OHMIC HEATING OF THE PROBE WIRE 


Probe wire heating due to the current through the wire is 
Qormic = IR, (38) 
where I is the current through the probe and R is the probe’s resistance. Brown in [Ref. 
14] measured this quantity to be 3.3x10°’ W at the small currents used in this type of 


probe, and we adopted his value for all of our modeling. 
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H. PROBE TEMPERATURE 


Ii 


Steady State Free Stream Velocity 


Solving the equilibnum equation for Tprote , using Eq. (15), the difference in probe 


temperature and atmospheric temperature AT versus altitude is shown in Figure 3-21. 
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Figure 3-21 : Differential probe temperature versus altitude 
with constant free stream velocity of 0 and 5 m/s with all 
probe areas constant at their maximum value for a 4.7 um plati- 
num coated tungsten probe. 


Normal Free Stream Velocity Fluctuations 


As the sonde rises, balloon pendulum motion and rotation occur. With the combined ver- 


tical horizontal and rotational velocity it is possible for the relative wind speed on the 


probe to drop to low values. The transient zero velocity condition is of interest in our 


model because forced convective heat transfer would approach zero. Near zero air veloc- 


ity, convection would be much smaller allowing the solar flux to dominate probe heating 


and increasing AT . In practice, though, free convection will create an induced effective 
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velocity over the probe wire, due to density differences near the wire so flow velocities 
will only approach zero. 
Figure 3-22 represents the effects of velocity fluctuations on the differential tem- 


perature as a function of altitude. 
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Figure 3-22: Simulated probe differential temperature for a 
normal velocity fluctuation of 2.5 m/s about a velocity of 5m/s 
versus altitude for the 4.7 um probe. 


3. Variation of Solar Illumination Area 


The effect on the differential probe temperature of varying the probe’s onentation 
with respect to the sun is shown in Figure 3-23 for sinusoidal change of area A2 as a func- 


tion of differential temperature and with constant altitude. 
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Figure 3-23: Effect of varying the solar illumination area, Ao, 
on the 4.7 um probe differential temperature with constant al- 
titude and air velocity (5 m/s). 


4. Probe Orientation with Respect to the Earth 


By sinusoidal variation of the area terms A, and As with periods of 5 and 20 sec- 
onds, the effect of relative Earth-probe orientations on the probe differential temperature 
for constant altitude and air velocity (5 m/s) may be calculated. Examples are shown in 


Figure 3-24. 
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Figure 3-24: Effect of sinusoidal vanation of Earth illuminated probe 
areas A, and As, (1.e., reflected solar radiation and Planck radiation 
from Earth) on the 4.7 um probe differential temperature with con- 
stant altitude and air velocity (5 m/s). 


As SUMMARY 


Relative energy contributions to the heat flux of the platinum coated tungsten wire 
are ranked by order of magnitude, and hence contribution to probe temperature, in Table 
3-3 for an air velocity of 5 m/s at sea level. Probe convection and direct solar radiation 
terms dominate, with convection being strongly dependent on free stream air velocity. 
Solar radiation effects are strongly dependent upon probe orientation. 

Fluctuations in the free stream air velocity over the probe can produce a maximum 
change in probe differential temperature of 0.1 K at 10 km. 

Direct solar radiation loading alone accounts for a maximum probe temperature 
change of approximately 0.111 K with constant free stream velocity, V=5 m/s at sea level. 


This change is a consequence only of a change in probe solar illumination area Ap. 
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Variation of the probe cross sectional areas presented to Earth (Ay and As) results 
in a maximum of 0.028 K variation of probe temperature with constant free stream veloc- 
ity of V=5 m/s at 30 km. 

Table 3-4 summarizes the maximum and minimum variation of the probe differen- 
tial temperature as a function of free stream air velocity, probe solar illumination area, and 


probe Earth illuminated areas. 


Direct Solar Illumination —— aa 
ARNOSONENE Planck Radiation en 


Table 3-3: Typical values for energy transfer mechanisms involved in the AFRL 4.7um 
platinum coated tungsten probe wire at constant air velocity of 5 m/s at sea level with 
maximum wire areas illuminated. 
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Table 3-4: Maximum and Minimum change in values for the AFRL 4.7 um probe differ- 
ential temperature for given parameter variation. 
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IV. EXPERIMENTAL MEASUREMENTS 


A. PROBE TEMPERATURE SENSITIVITY CHARACTERISTICS 


1. Purpose 


Verification of the actual probe response charactenstics when subjected to solar 
radiation and varying wind speeds was done experimentally to validate probe modeling of 
the solar and convective heat transfer mechanisms. In this experiment we measured the 
differential temperature of the probe in response to one solar flux of radiation as a function 
of wind speed. The effect of Earth reflected sunlight at zero wind speed on the probe was 


also approximated experimentally. 


2: Experimental Procedure 


The platinum coated tungsten probe was mounted in a wind tunnel with air flow 
perpendicular to the wire (Figure 4-1). A Radiac ELH 300 W tungsten halogen lamp was 
focused uniformly on the probe through an optical lens, and the lamp power was adjusted 
to correspond to one solar flux after being chopped with a Laser Precision Corporation 
CTX-534 Radiation Chopper. Calibration of our unit solar flux source was done by direct 
comparison with the Sun at local solar noon on a clear day at latitude 40N, as measured 
by a United Detector Technologies model 255 silicon photodetector. A uniform spectral 
response measurement was made with a Scientec 362 Thermopile Power Energy meter to 
cross check the spectral response of the photodetector. To the accuracy of our measure- 
ments (about 1%), the halogen lamp flux and the solar spectral flux produced the same 
probe heating with both the silicon and broadband detectors. The signal was passed 
through a Stanford Research Systems SR560 Low-Noise Preamplifier and voltage read- 
ings were then taken using a Hewlett-Packard 3561A Dynamic Signal Analyzer and Ni- 


colet Pro 30 digital oscilloscope for wind speeds from 0 to 6.0 m/s. Wind velocities were 
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measured using an Omega HHF 52 hot wire anemometer accurate to 0.1 m/s. Three sepa- 
rate experiments were conducted on 24-25 Apmil, 16 May, and 29 May 1997 using three 
different platinum coated tungsten probes. Probes #1 and #2 were from similar batches of 
platinum coated tungsten wire while Probe #3 was a much older (1989) probe from a 
separate production run of wire obtained from the Air Force Geophysics Directorate. 

The effect of Earth reflected sunlight was measured using a 41% reflective gray 
surface behind the probe. The reflected sunlight measurements were made at zero wind 


speed. 


3. Results 


Experimental results are tabulated in Tables 4-1 to 4-5. The rms voltages read 


from the HP-3561A spectrum analyzer were converted to peak-to-peak voltages by multi- 


plying by a factor of 2/2 and dividing by the gain applied at the SR560 preamplifier, 
since the wave form of the chopped light was nearly sinusoidal. For the second experi- 
ment (16 May 97) this procedure was modified slightly to include a correction factor, of 
2.22, to account for an observed non-square, chopped wave form from the halogen lamp 
source. The probe differential temperature was obtained by solving for AT in the equation 
for probe resistance [Ref. 13:p. 11] 

R, = R, (1+ aAT), (39) 
where Ro is the cold reference probe resistance. A different probe was used in each ex- 
periment so this value was measured as 19.94 Q (24-25 April 97), 18.81 Q (17 May 97), 
and 25.68 02 (29 May 97). The constant alpha was assumed to be a=0.0042 for platinum 
coated tungsten wire [Ref. 30]. 

An upward drift in measured voltages was noted after continual illumination of 
about 60 s or more. We suspect that this was most probably due to probe support heating 
effects as shown in Ref. 5:p. 16. Severe direct thermal contamination was noted with the 
halogen lamp closer than about 0.5 m, obviating the use of a lens to separate the hot lamp 


from the wind tunnel. 
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The observed temperature change depended on the chopping frequency and the 
probe response time. The probe response time varied between 0.6 to 2 ms and acted like 
a one pole low-pass filter at these frequencies. Collecting measurements at low frequen- 
cies, on the order of 5 Hz, was not always possible because of interference from ambient 
turbulence. In these cases higher frequencies around 25 Hz were used after including a 
correction factor. 


Probe #3 showed a lower probe differential temperature (Table 4-4) than Probes 


















#1 and #2. Possible explanations for this discrepancy are investigated in Chapter V. 
Differential Tem- 
perature (K) 


Air Velocity Measured Volt- Peak-to-Peak | Actual Volt- AR (mQ) 
(m/s) age hp 3561A Voltage - age (uV) 
mV- rms Corrected 
129.0 364.9 0.095 [0.120] 
Oma 1143 | 323.3 0.085 _[0.091 


m 
Cm . 
0 
0.071 
0.061 
0 . ; 
0 





Le 

Zr 

a 
0.059 
0.055 
(60 [67.001 189.5 3.79 4.14 0.049 
Table 4-1: Experimental probe differential temperature measurements of Probe #1 (new) 
for one solar flux from a halogen lamp using a chopping frequency of 50 Hz and cor- 
rected temperatures for 0 and 1 m/s air velocity measurements (19 April 1997). 


Measured Volt- Peak-to-Peak Differential Tem- 

| (m/s) age hp 3561A Voltage - age (pV) perature (K) 
mV- rms m Corrected 

(Oo 1281 I 3 a Eee 

5.0 | 

8 

5 4.6 




















le 

on 

3. , 

4. 5.0 10.060 

5 76.8 | 217.3 
46 10.054 
Table 4-2: -Experimental probe differential temperature measurements of Probe #1 (new) 
for one solar flux from a halogen lamp using a chopping frequency of 25 Hz and cor- 
rected temperatures for 0 and 1 m/s air velocity measurements (24 April 1997). 
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m/s Pro 30 (mV- rms age ture (K) 
a Pee eer ee 030 ae 


Table 4-3: Expernmental probe differential temperature measurements of Probe #2 (new) 
for one solar flux from a halogen lamp using a chopping frequency of 5 Hz (16 May 
1997). 


m/s 3091 (mV- rms age ture (K) 
2 See h 


Table 4-4: Experimental probe differential temperature measurements for Probe #3 
(1991-3.8 um) one solar flux from a halogen lamp using a chopping frequency of 10 Hz 


(29 May 1997). 
Air Velocity (m/s) | 25 Hz, Measured 5 Hz, Measured 
Voltage (m Voltage (m 


Table 4-5: Correction factor for low air velocity and higher frequency chopping rates. 


















Adding a gray cardboard background for reflection of the halogen lamp provided 
the probe with a 27 steradian solid angle reflective background surface which roughly 
simulated that of the Earth. Table 4-6 shows the increased signal differences due to 
“Earth reflected solar radiation” from the absence of a reflective surface, gray cardboard, 


reference (white) surface, and a first surface mirror. 
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Percent Change (%) 














3561A (mV- rms 
9 
[White (00%) ——~=«dt6 
| First Surface Mirror (~100%) |60.8 = = «|. 43 
Table 4-6: Experimental data for the effect of reflective backgrounds on differential probe 
temperature of Probe #2 (new) with one solar flux of direct radiation at zero air velocity, 
sea level, and chopping frequency of 5 Hz (29 May 1997). 






> | 
N 





Sources of error for this experiment include the signal lost due to the high chop- 
ping rate in the earlier data. This is partly resolved with the correction factor in Table 4-5. 
Other minor sources of error involve probe support heating effects, inaccuracy in wind 
velocity measurements, and assumptions that the halogen lamp accurately reflects the 


spectral output of the sun. 


4, Summary 


At sea level, with no forced air flow over the probe, radiation heating countered by 
free convection and Planck radiation caused probe heating about 0.130 K above ambient 
temperature. This differential temperature decreased as forced air flow over the probe 
effected more efficient forced convection cooling to about 0.055 K at 5.0 m/s. 

Earth reflected solar radiation contributed an additional 19% increase in probe 
temperature at zero forced air velocity and would be expected to remain independent of 


air velocity and approximately independent of altitude. 
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B. MEASUREMENT OF PROBE WIRE RESPONSE TO DIRECT SOLAR 
RADIATION 


a’, Purpose 


The probe response to direct solar illumination at sea level was measured under 
natural conditions to detect any additional error sources or procedural omissions in our 


laboratory measurements and simulations. 


ze Experimental Procedure 


Using a similar procedure to that described in Section IV-A, above, the platinum 
coated tungsten probe was subjected to direct solar illumination on the roof of the NPS 
physics building (Spanagel Hall) at approximately solar noon on a clear day. High back- 
ground turbulence forced the use of a spectrum analyzer with chopping frequencies of 25 
Hz (see Figure 4-2) but otherwise the experimental procedure followed that of Section 
IV-A. 


3. Results 


The results of this experiment are tabulated in Table 4-7. Multiple data points for 
each air speed where taken and averaged because of the difficult measurement conditions. 

Possible error sources include effects from heated air rising from the rooftop con- 
crete surface. This caused noticeable fluctuations as the fan pulled the hotter turbulent air 
over the probe. The wind speed on the roof varied from zero to almost 5 m/s. Although 
the experiment was sheltered from direct wind, some breeze was evident. Ambient at- 
mospheric turbulence made it difficult to collect probe solar illumination data accurately. 
A chopping frequency of 25 Hz facilitated in this but forced a correction to the results as 


discussed in the previous section. 
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4. Summary 


Probe temperature rises, as seen in Table 4-7, from direct solar illumination were 
about 30% lower than those for the laboratory measured halogen illumination, after ap- 
plying the correction because of the higher chopping frequency. The error was most likely 
introduced from ambient air flow fluctuations encountered on the rooftop. Considering 


these corrections and the unsteady air, agreement with experiment was good. 












Actual | Rw (mQ) AT (K) 


Voltage 






Air Velocity | Measured 
Voltage 


Average AT 





| [0.0980] 
0.0749 


0.0462 
0.0430 


0.0337 


Table 4-7: Data from rooftop experiment using a chopping frequency of 25 Hz (16 May 
97). 


GC: MEASUREMENT OF PROBE WIRE VIBRATION EFFECTS 


1 Purpose 


Historically, error sources in the balloon thermosonde data were thought to result 
partly from probe wire vibration effects. As air flow over the probe induces vibration in 
the probe fine wire filament, it would be subjected to periodically varying solar radiation. 
This vibration, along with the package pendulum motion and rotation, could induce a sig- 


nificant artifact in the data. 


2. Experimental Procedure 


With the probe mounted in the wind tunnel (Figure 4-3) the behavior of the thin 
wire filament was investigated with respect to air flow speed and orientational depend- 
ence on vibration. Illumination of the probe wire was by a Hughes 7.35mW helium-neon 
laser beam chopped at 90 Hz. The probe output signal was amplified by a Stanford Re- 
search Systems SR560 Low-Noise Preamplifier, passed through a Wavetek Brickwall 
bandpass filter (model 753A), and read on both a Hewlett-Packard 3561A Dynamic Signal 
Analyzer and a Tektronix 2236 100 MHz oscilloscope. The bandpass filter was set at 80- 
100 Hz and the gain was optimized for the signal being read. The probe was also visually 
monitored using an American Optical microscope placed above the probe but out of the 
air flow. 

Turbulent flow over the wire was induced using a thin rod placed just to one side 
and upstream of the probe. Reynolds number for the flow over this rod is on the order of 
Re=1000, which is in the turbulent regime and corresponds to a vortex shedding frequency 
of 50 Hz from [Ref. 31:p. 15] 


SV 
Te (40) 


with S the Strouhal number ( S=0.15 in this case), V the free steam velocity, and D the 


thin rod diameter. 
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3. Results 


Models for flow induced vibration of thin wires developed by Perry in [Ref. 13:p. 
99] and Blevins [Ref. 31:p. 15] predict vibration frequencies about two orders of magni- 
tude greater than what were observed visually. Because of the very low Reynolds num- 
bers involved with the platinum coated tungsten wire, the probe wire vibration is most 
likely due to turbulence in the free stream flow over the wire and not to resonance or 
vortex shedding effects from the wire itself. As previously shown, the Reynolds number 
throughout the realm of our model is less than 4 (Figure 3-6 and 3-7), which is under that 
for separation of flow and turbulence for a circular cylinder which occurs with Re>S [Ref. 
21:p. 52]. Observed vibration of the wire was very dependent on the presence of turbu- 
lence in the air flow (necessitating the thin rod obstruction) and also showed a marked vi- 
bration onset near velocities of 2 m/s with oscillations continuing at nearly constant 
frequency and magnitude to the limit of our test, 5.8 m/s. 

With the probe oriented such that the wire was vibrating in and out of laser illumi- 
nation, the presence of filament vibration had a noticeable effect on the output signal char- 
acteristics as observed on the oscilloscope. Fluctuation about the mean sinusoidal 
(chopped) output would increase with visible turbulence but was sufficiently random and 
not of a magnitude that would allow spectral components or harmonics to be noted above 


the background noise on the Hewlett-Packard signal analyzer. 
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4. Summary 


Although the probe wire was noticeably vibrating in and out of the laser beam at a 
frequency of less than 50 Hz, near that of the thin rod vortex shedding frequency, no ap- 
preciable signal could be quantified to the accuracy of our experiment. From this we con- 
clude that probe solar illumination vibration effects are not a significant source of data 
contamination, although they can occur in the presence of velocity fluctuations in the air 
flowing over the probe. On the other hand, airflow over the probe will displace the probe 
wire, since it is not taught. Differential displacement of the probe wires, in a dual probe 
configuration, will alter the area exposed to the sun. The magnitude of this displacement 
depends on the drag force, which declines with altitude. This reduction with altitude can 
produce the pattern seen in Figure 2-5, where the solar induced contamination declines 


with altitude above 20 km. 


D. MEASUREMENT OF COEFFICIENT OF THERMAL CONDUCTIVITY 


1. Purpose 


Four analytical methods based on experimental data for the modeling of the coeffi- 
cient of thermal conductivity, h, have been presented in Chapter II]. To provide a check 
of the analytical expressions we chose to independently measure h for our probe. Two 
experimental methods were used. Atmospheric pressure was varied from sea level to 1 
torr, with air velocities over the probe ranging from 0 to 5.8 m/s. Three different probes 
were measured with tungsten coated platinum wires from two different manufacturing 
batches. Probes #1 and #2 were purchased new for this work while Probe #3 was previ- 


ously used in 1991. 
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Z, Experimental Procedure 


The platinum coated tungsten probes’ response to a square wave input was meas- 
ured on both a Nicolet Pro 30 and Nicolet 3091 digital oscilloscope (Figure 4-4). Our 
first method used the heat rate of change equation given in Richards [Ref. 22:p. 87] 


or OT 
IP = I(T ype — T) + kB ~ 208, (T* ~ Te), (41) 
where r is the probe radius, and To is the effective black body temperature of the sur- 
roundings. Neglecting the last two terms in Eq. (42) (conduction and radiation terms) the 
estimated rise time of the probe response is: 


C 
4 pDe 


°= GDI) 4h” | " @ 





where T is the rise time, c, is the heat capacity (133 J/kg-K), D is the probe wire diameter 
(m), and p is the density of tungsten (19300 kg/m’) [Ref. 20,p. 543]. 
The second technique relied on measurement of the asymptotic temperature rise 


from Joule heating of the probe wire and the relations: 





V? 3 
q=— = R=RAAT, (43) 
R(T) = R,(1+ @AT), and (44) 
>= VLR, 45 
probe — V, -V, : ( ) 


Here V;, is the voltage applied by the function generator with series resistance R;, V2 1s the 


measured voltage output on the oscilloscope, and Rpyote is the probe resistance (Figure 4- 


5). 
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Figure 4-5: Circuitry for the measurement of probe resistance. 


Air velocity was varied from 0 to 5.8 m/s for both of the above methods and the 
minimum voltage possible was used to reduce overheat effects. Conduction effects 
through the probe wire supports were measured as described in Lomas [Ref. 19:p. 69] by 
placing the probe in a vacuum and measuring the heat transfer rate to the probe support 
needles by: 

F move = IR. (46) 
Finally, one probe was subjected to a vacuum under a zero air flow condition and varying 


pressure as measured with a Vaisala PA 11 digital barometer. 


3. Results 


The three sets of data are shown in Tables 4-8 to 13. Values for h were calculated 


using the MATLAB program listed in Appendix D. 
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Sensitivity to probe overheating was found and is shown in Tables 4-13 and 4-14. 
This overheat condition was thus minimized as much as was possible when making meas- 
urements. The reason for this sensitivity is unknown and is still being investigated. 

Contamination of the probe wire through normal use resulted in an increase in the 
observed probe response time by about 20% over that of a thoroughly cleaned probe. 
This effect was thought to be negligible and was first discounted as an error source. Perry 
in [Ref. 13:p. 56] discusses its possible causes and prevention. 

Conduction of heat through the probe wire supports was measured to be 0.0015W 
at an overheat of 53 K, zero wind speed, and a pressure of 0 mBar. This power flow de- 
creases as the square of the input voltage and is therefore very small for our voltages and 
corresponds to a 2.8x10° W/K temperature difference between the probe wire and the 


probe supports. 


_—— (mV) al (ms) 
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Table 4-8: Experimental data for time constant measurement technique of Probe #1 (new) 
using Nicolet 3091 and varying air velocity(1 May 97). 
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Table 4-9: Experimental Data for time constant measurement technique of Probe #2 (new) 
using Nicolet 3091 and varying air velocity (14 May 97). 
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Table 4-10: Experimental Data for time constant measurement technique of Probe #3 
(1991) using Nicolet 3091 and varying air velocity (29 May 97). 
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Table 4-11: Experimental Data for asymptotic temperature change technique of Probe #1 
(new) using Nicolet 3091 and varying air velocity (2 May 97). 
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Table 4-12: Experimental Data for asymptotic temperature change technique of Probe #3 
(1991) using Nicolet 3091 and varying air velocity (29 May 97). 
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Table 4-13: Experimental Data for asymptotic temperature nse technique for Probe #1 
(new) using Nicolet Pro 30 and varying pressure and zero wind speed (5 May 97). 


Input Volt- 1/e Rise Time | Overheat (K ) | h (W/m?-K) 
age (V p-p ms 


279.2 286.8 0.750 6475 
2508 | 32 OO 20 8120 


Table 4-14: Sensitivity of coefficient of heat transfer to overheat condition (8 May 97). 





















4. Summary 


The values for the coefficient of heat transfer, h, using the mse time method are 
consistent with those associated with the calculated ATs from Section IV-A. The values 
of h derived from the time constant method, Tables 4-8 and 4-10, are almost one half 
those from the asymptotic temperature rise technique, Tables 4-11 and 4-12. 

The values of h are very sensitive to the amount of overheat, and therefore cur- 
rents through the probe should be minimized to reduce error. Although measurements 
were made with four significant figures it was hard to achieve 20% accuracy in the values 


of h. 
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E. MEASUREMENT OF PROBE WIRE RESISTANCE 


1. Purpose 


The previous section measures resistance changes with probe current dynamically 
using digital oscilloscopes. Another method of determining h measured the probe resis- 
tance change with varying input current and changing wind velocity for static changes in 
the resistance. This version of the method in Section V-D relies on the accurate input of 
current to the probe and precise determination of current through the circuit using meas- 


ured voltage across the load resistor, Figure 4-5. 


De Experimental Procedure 


As shown in Figure 4-6 the probe resistance was measured using a Fluke 8840A 
multimeter and a 5’/2 digit Keithley 197 Autoranging Microvolt DMM _ in response to a 
current input from the Hewlett Packard 6235A power supply. Different supply current 
levels were tried in order to minimize the overheat effects and change in load resistance 
with temperature. Relations for the calculation of AT and the coefficient of heat transfer 
have been developed in the previous section. Contact resistance in the probe was meas- 
ured and considered in the results. Accurate determination of the load resistance was 
done with three independent instruments using the four wire measurement technique [Ref. 


36]. 
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3. Results 
The results using 
P 2k ae (47) 
are shown in Table 4-15. The contact resistance of the probe was measured as 0.488Q 
and has been subtracted from the tabular value of R2 shown. 
Possible error sources include a thermocouple type response of the probe wire 


supports and some heating effects from the Hewlett-Packard 331210A function generator. 


Air Ve- | V; (Volts) | R; (Q) I, (mA) | V; (Volts) R; (Q) AT (K) h (W/m’K) 
| locity | with 2 
(m/s) Ro=19.285 
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Table 4-15: Results of probe resistance method for determination of the coefficient of heat 
transfer of Probe #2 (9 May 1997). 













4. Summary 


Values for this method, which is a refinement of that used for Table 4-11 and 4-12, 


are also higher than predicted from either our model or the rise time technique. 
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Mi DISCUSSION 


A. MEASURMENT OF PROBE WIRE DIAMETER AND COMPOSITION 


1. Purpose 


All of the preceding experimental measurements and simulations rely on, among 
other parameters, the accurate knowledge of the platinum coated tungsten probe’s di- 
ameter. The manufacturer, Thermo Systems Incorporated (TSD, gives this as 3.8 um, but 
data taken in the course of this research raised doubt as to accuracy of this number. 

Also questioned during the course of our research were the surface characteristics 
of the platinum coated tungsten wire as they relate to its reflectivity and heat transfer 


characteristics. 


2s Experimental Procedure 


A new sample purchased directly from TSI in 1997, and an older one obtained 
from the AFRL in 1991, were placed in the chamber of a Cambridge 20 KeV scanning 
electron microscope (SEM). The SEM had recently been calibrated and provided a scale 
accuracy of about 10 nm for our experiment. Measurements of the diameters were then 
taken at several positions along the wires, along with images of the regions measured. 
Images of the probe wires were also stored to disk for later analysis. A second sample of 
a Dantec pure platinum 5 um wire was also imaged to provide further information about 
other vendor wire surface characteristics. 

A Kevex Energy Dispersive X-ray (EDX) spectrograph then took representative 
measurements of the elemental composition at locations of the newer wire. These were 
recorded both where the platinum coating was present and at locations where spalling, 
thermal cycling, and possibly other processes had removed large patches revealing the 


underlying bare tungsten. 
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Sh Results 


Conservative measurements of the wire diameter at several locations, Figures 5-1 
and 5-2, for the new wire revealed an actual diameter of 4.7 um. Measurements of the 
older wire sample Figure 5- 3 showed its diameter to be 3.7 um. 

The characteristics of the new wire revealed an average of up to 22% bare wire 
(absence of platinum coating) over most of the surface with several large patches of bare 
tungsten, Figure 5-4 and 5-5. Striations in the platinum coating, very similar to the bark 
on aredwood tree, created a very rough and non-uniform reflective and radiative surface. 
(Figures 5-1,2,4,5,and 6). 

In contrast to the rough surface seen in the new wire, Figures 5-1] and 5-2, the 
older wire’s surface is significantly smoother and more reflective, Figure 5-3. Although 
the image quality is not quite as good as those for the newer wire, the lack of deep stria- 
tions and dark patches is still discernible. The much smoother surface texture of a Dantec 
P11 pure platinum probe is shown in Figure 5-7. 

Figure 5-8 shows that the tungsten areas were relatively pure in composition, while 
Figure 5-9 shows the significant amounts of copper that were detected by the EDX spec- 
trum analyzer in the platinum areas. This probably indicates the presence of either con- 
tamination or impurities in the wire material, possibly introduced during the 


manufacturing process. 
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Figure 5-5: SEM image of striations and small bare patches (30 May 97). 
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Figure 5-6: SEM image of wire end. 
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Figure 5-8: X-ray emission spectral components of bare tungsten wire. 
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Figure 5-9: X-ray emission spectral components of platinum wire coating. 
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4. Summary 


The new sample of the TSI platinum coated tungsten probe wire diameter was 
measured to be 24% greater than quoted by its manufacturer, while the older wire ob- 
tained from AFRL was 3% smaller than stated. Surface characteristics of the new wire 
showed a darker surface than appropriate for pure platinum and thus will introduce an er- 
ror in the values for the emissivity and albedo assumed in the solar heating model calcula- 
tions. Other heat transfer related effects from the rough and flaked surface of the new 
wire batch cannot be determined within the scope of this thesis research. 

Overall, the older platinum coated tungsten wire showed very similar surface char- 
acteristics to the pure platinum wire, as expected. The newer wire was significantly differ- 
ent and showed marked surface roughness, bare patches, striations, and significantly larger 


diameter than specified. 


B. COMPARISON OF MEASURED AND MODELED COEFFICIENTS OF 
HEAT TRANSFER, h 


Figures 5-10 and 5-11 compare the coefficient of heat transfer as modeled in the 
MATLAB program of Chapter III, for Probes #1 and #3 and the two different measure- 
ment techniques, i.e., time constant and asymptotic temperature change. The large dis- 
crepancy using the asymptotic temperature change technique and the time constant 


procedure for both probe wires is not understood and is still under investigation. 
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Figure 5-10: Modeled and measured h, for the 3.8 ium Probe #3 for 
measured rise time and asymptotic temperature increase techniques. 
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Figure 5-11: Modeled and measured h, for the 4.7 um Probe #1 for meas- 
ured rise time and asymptotic temperature increase techniques. 
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C; COMPARISON OF MEASURED AND MODELED DIFFERENTIAL 
TEMPERATURES 


Figures 5-12 and 13 compare the differential temperature rise from solar and oh- 
mic heating effects, measured experimentally in Section IV-A, to temperatures computed 
using the MATLAB model of Chapter III. Measured Probe #3 data , in Figure 5-12, 
compares favorably to the models for emissivities of 0.435 (platinum coated tungsten) and 
0.444 (pure tungsten). The 4.7 um Probe #2, Figure 5-13, does not correlate as well most 
likely due to the unknown heat transfer and absorption effects caused by the poor surface 


characteristics previously mentioned. 


----R38034 | 
|—— Prebe#B 


E ----RW380.43 


Air Velocity (ms) 





Figure 5-12: Modeled (dashed and dotted)and measured (solid) data for the 3.8 um plati- 
num coated tungsten wire Probe #3 differential temperature from solar and ohmic heating. 
Emissivities of 0.334 (pure Pt), 0.435 (Pt-W), and 0.444 (pure W) are shown for the 
modeled curves. 
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Figure 5-13: Modeled (dashed and dotted) and measured (solid) data for the 4.7 um 
platinum coated tungsten Probe #2 differential temperature from solar and ohmic heating. 
Emissivities of 0.435 (Pt-W) and 0.444 (pure W) are shown for the modeled curves. 


D. EFFECT OF PROBE WIRE DIAMETER VARIATION 


The MATLAB simulation program developed for Chapter III provides a way to 
investigate the effect of different probe wire diameters on the differential temperature. 
Changes in the probe wire diameter with altitude on the differential temperature are re- 
flected in Figure 5-14. Simply reducing the wire diameter will reduce the probe’s differ- 
ential temperature. For example, a reduction in wire diameter from 5 to 2 um would 


result in a AT reduction from 0.16 to 0.07 K or a 64% improvement, at 30 km. Further 
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reduction of probe diameter below 2 um into the more pronounced non-linear region 


would have a larger effect. 


0.35 


oO 
w 


© 
NO 
On 


0.2 


0.15 


Differential Probe Temperature (K) 


0.05 





Probe Wire Diameter (microns) 


Figure 5-14: Effect of varying the AFRL probe wire diameter for 
the platinum coated tungsten wire on differential temperature at con- 
stant altitude and air velocity of 5 m/s. 


E. EFFECT OF PROBE MATERIAL AND DIAMETER CHANGES 


In light of the very poor surface characteristics discovered for the TSI wire, other 
probe wire materials were considered. The dominant parameter affecting the probe wire’s 
solar absorption is the spectral emissivity. The spectral emissivity of several favorable re- 


placement materials where calculated using the following relation: 
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éy - co ? (48) 
| pay 
0 
where a(v) is the absorptivity as a function of frequency and 
8zv° hv 
Av\dv = 5-5, —av, (49) 
e™ —] 


which is the Planck blackbody spectrum [Ref. 32:p. 17]. Curves for the spectral emissivity 
versus energy from values tabulated in the CRC Handbook are compared with the corre- 
sponding black body spectrum of the sun, a 6000 K body, in Figure 5-15. The black body 
spectrum was used as the weighting factor for Figure 5-16, which compares the integrand 
of Eq. (48) for the four metals. Completing the above integrals, over the range corrte- 
sponding to the transparent region of the atmosphere, about 0 to 195 nm, we then arrive 
at Table 5-1 which lists the spectral emissivity for the five wire materials investigated. 
Revisiting the previous section and using the spectral emissivities of the new mate- 
nals we can plot the change in differential probe temperature with varying wire diameter at 


a constant altitude and air velocity of 5 m/s as shown in Figure 5-17. 
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Figure 5-15: Emissivity versus energy for some probe matenials of interest 
with 6000 K black body for comparison. 
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Figure 5-16: Weighted emissivity curve for different probe wire matenials. 
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Figure 5-17: Effect of change in spectral emissivity on differential 
temperature and wire diameter at 30 km and 5 m/s. 
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Spectral Emis- 
0.444 0.334 0.116 0.073 


Sivity 
Table 5-1: Average spectral emissivities of some selected metals for the solar spectrum. 









In Figure 5-17 it is evident that a change of wire material from the present plati- 
num coated tungsten to silver would improve wire solar absorption performance by 67% 
with both wire diameters at 5 um and at an altitude of 30 km. Reduction of the wire di- 
ameter to 2 um and change in material to silver would reduce the solar heating by 85% at 
30 km. Since Cy’ is proportional to AT’ a factor of 6.7 reduction in the solar induced 


temperature would reduce Cry’ error by a factor of 45. 
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VI. CONCLUSIONS 


In this thesis, the measured values for the solar heating of the 4.7 um platinum 
coated tungsten probe varied from 0.055 K at 5 m/s to 0.12 K at 0 m/s and 0.034 K at 5 
m/s to 0.085 K at 0 m/s for the 3.8 um probe, at sea level. These values were 12% to 
34% larger than predicted by the detailed model for the 4.7 um probe. Much better 
agreement occurred for the 3.8 um platinum coated tungsten probe. Measurements and 
model results agreed to within 1-7%. Direct solar heating of the probe is significant. 
Temperature fluctuations over 0.175 K, at 30 km, can be induced for both 3.8 and 4.7 um 
probes through the periodic modulation of the rotation and pendulum motion of the sonde 
package both due to the direct illumination (heating) of the probe but also because of the 
changing air flow (cooling) over the probe (Table 3-4). Notably, the maximum change in 
the probe’s differential temperature occurs at an intermediate altitude of 11 km with de- 
creasing differential magnitudes occurring above and below this altitude. This effect is 
due to solar heating and is best shown in Figures 3-21 and 3-22. 

The uncertainty in measured versus modeled values for the differential temperature 
of the 4.7 um platinum coated tungsten probe in the atmosphere is most likely due to the 
complicated surface characteristics of the TSI wire. Convection, radiation, and air flow in 
and around this rough surface introduces unknowns which could easily be of sufficient 
magnitude to account for the measured differences. 

The determination of the value of the coefficient of heat transfer, h, is important in 
characterization of the probe’s ability to transfer energy and was therefore modeled and 
measured with great care. Final values of h calculated through the model and measured 
experimentally agreed to within 20% at air velocities of 0 m/s and 10% at air velocities 
near 5 m/s. Probe wire surface characteristics most likely account for the consistently 
higher differential temperatures, of about 0.002 K, measured experimentally for the 4.7 


uum probe. 


wo 


Susceptibility to the direct solar heating effects can be substantially reduced by 
minimizing the probe wire diameter. Additional improvement to solar contamination ef- 


fects can be realized by changing the probe matenial, notably to silver or aluminum. 
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VU. RECOMMENDATIONS 


The best solution to solar contamination of data would be to perform only night 
data collection. As this is not always an option in all research and tactical situations other 
options are necessary. The TSI platinum coated tungsten wire measured has been shown 
to have very poor characteristics, from a solar absorption standpoint, and thus would be a 
prime candidate for replacement. This may have been a result of manufacturing defects. 
Options include pure platinum wires which present a much more reflective and thermally 
smooth surface or another material such as silver or aluminum. The diameter of the probe 
wire is also a major contribution to solar contamination and should be reduced. Similarly, 
the thermal heating of the probe supports must be minimized using highly reflective mate- 
rials. 

Replacement altogether of the probe assembly has been proposed. Some replace- 
ments include an acoustic sounder sensor or a laser/optical device. If the old design probe 
must be used for balloon sounding, further research could be performed into stabilization 
techniques for the sonde to minimize or remove the geometry dependence of the balloon - 
sun system. For the near term, reduction in the probe wire diameter and shift to a new 


wire material, such as silver or aluminum, is highly recommended. 
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APPENDIX A- AIRBORNE LASER PROGRAM SUMMARY 


Table A-1: Highlights of the Airborne Laser Program acquisition process [Ref. 6:p. 12). 
September 1992 Boeing and Lockheed receive separate contracts to assess how well existing large air- 
planes, the 707, 747, 767 and B-52 would perform while carrying a high-energy laser 
and beam control system. Both teams separately reach the same conclusion and rec- 
ommend a 747 platform to the Air Force as the best option for the mission. 







| 1993 Team ABL ————— s+). | submitted 


_— 1994 (MSO) Team ABL eee Rockwell selected as two competing teams for Concept Study Design 
($21 million, 33-month contract) 


March-November 1995 Boeing conducts series of “Flex Blue” tests to determine vibration environment for the 
Airborne Laser. 


Apnil 1995 TRW begins construction of BDL-2, (building block COIL). 
May-October 1995 Boeing Wind Tunnel Testing of Turret Designs. 
July 1995 TRW TRACE facility produces first laser light to demonstrate 10 kW laser power goal 


—— 1995 Air Force Operational ee Review (OCR). 
October as Lockheed Martin beam control An ellen + pointing and adaptive optics compensation 
demonstrations 
| anuary 1996 Draft RFP for PDRR _————_ ... released 


February 1996 Lockheed Martin beam control system brassboard operational 


Saar Sar TS 


November 12, 1996 (MSI) | Air Force selects Team ABL as one contractor for Airborne Laser Dem/Val. $1.1 billion 
5 year Program Definition and Risk Reduction (PDRR contract to Team ABL 
(Boeing/Lockheed Martin/TRW) (cl) 


1996-2000 PDRR Phase I (build one PDRR aircraft). 


2000-2001/2 PDRR Phase II (test and fly PDRR aircraft. 


2001/2-2004 (MSI) Engineering Manufacturing and Development (EMD). 
2004-2008 (MSI) Production (IOC 3 aircraft 2005/6, FOC 7 aircraft 2007/8). 
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Figure A-2: Airborne Laser Program acquisition timeline. 
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Figure A-2: ABL Program structure 
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APPENDIX B - PROBE SIMULATION PROGRAM 


$Daniel J. Richardson 
$13 May 97 


&%Probe Temperature vs. Air Temperature 
Using Free and Forced Regimes 
$and Low Re number Equations 


%Set Constants: 
alpha=0.435; 
Gpsi1con—0.065; 
gamma=0.39; 
Sigma= 5.67e-8; 


$Solar Absorptivity of Probe 

$Emmisivity for 300K Radiation 

$Solar Energy Reflected by the Earth (Albedo) 
&Stefan-Boltzman Constant [W/m*2.k%4] 


Pe=615; Long Wave Radiation from Earth [W/m%*2] 

Pa=250; Long Wave Radiation Emmited by the Atmosphere [W/m%2] 
k=0- 3; Atmospheric Transmission Term 

D=4.7e-6; Diameter of Probe Wire [m] 

L=4.0e-3; Length of Probe Wire jm] 

v=5; $Balloon Ascent rate 

betal=1.458e-6; Dynamic Viscosity constant [(kg/smk%*1/2] 

S=1 10. 4: $Sutheriand's Constant [Kj 

zO0=8.1; $Scale Height [km] 

pO=101300; $Scale Pressure [Pa} 

M=0.0289644; Mean Molecular Weight [kg/mot]} toa? 
R=8.31; Universal Gas constant [J/mole K)} 

Psurf=890; Solar flux at the surface of the Earth [W/m*2] 
W=2.5; Wind Velocity Factor 

c=1010; $Specific Heat [J/kg K] 

G=9-.62; $Gravity 


beta2=3.66e3; 


$Coefficient oc thermal expansion [1/K] 


Atmospheric Variables: 
Altitude=0:10:30e3; 
Kt=Altitude; 
Ps=Altitude; 
A2=Altitude; 
t=Altitude; 
Pr=Altitude; 
Re=Altitude; 


$Altitude in {m} 

SDefine array size 
$Define array size 
$Define array size 
$Define array size 


for i=l:length (Altitude); 
emer ol cue (i) <1 .00G) 
TA y=(Albereude (i t71000=44762) 7 (-0, 15443); 
Ket(iy=( (Altitude(1)/1000)—S02 57 )7 (=19e3247); 


$Temperature of Atm. [K]} (O-l1lkm) 
$Thermal conduct. [W/m-K] (0-12km) 


elseif (Altitude(i}>=11000 & Altitude(i)<20000) 


Mer r= 21.2 % és (11-20km) 
Ken i) =O. O95: % " (12-20km) 
else 
Vi) —(Aleautcuae( i )/1000+142:.75)/0.7S: % ie (20-30km) 
Kea ta ercude (a1) /1000+T206.5)/1 1636; % i (20-30km) 
end 


end 


&Solar Flux [W/m%*2] 
$Pressure [Pa] 

Dynamic Viscosity 

$Density 

Kinematic viscosity [m%2/s] 


Ps=(Altitude.°0.584)+890; 
P=p0*exp( (-Altitude/1000)/z0); 
mu=bpetal*(T. 7 (37 2) ) 2/4T+s) ; 
rho=(M.* F) o/s a 
nu=mu./rcho; 

$V=W* randn(size(Altitude) )+5; $Random Air velocity generation around probe 
deltaT=0.1; 
Gr=('(rho.°2)*Gq*betaZ*deltalt™ (b-3))../mu; 
Re=D*abs(v)./nu; 
Reefr=(Re.726 (0. 92Gr. 6.418) 472). 0. 5; 
Pe=(Mu.* 6) 27 Ke, 


$Grashof number 
$Reynolds Number 
tEffective Reynolds number mix.Convection 
Prandtl Number 
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for j=l: lengun{Alcicude, 


if (Reiq)<0-Cils2) Check for Free Convection ONLY regime 
ht} )=(0e32* kets Sp.131 Bayan 
else 
hy }={0.24+0.56* Reet (a7: (024 5eeeke (7) 27 0); 6  SMixed convection 
end 
ena 
Nu=D*h. 7 ke; Nusselt Number 


Areas Seen by the Probe: 


Al=pi* Dae. Surface Area of Probe wire 

A2=(D*L)- Area of Probe seen by Sun 

t=0:@225. 750; 

ps=randn(1); $ohase factor for Sun area 

pe=randn(1); phase factor for the Earth area 

Wp=(2*paniy 20; $angular frequency of Pendulum motion 

Wr=(Z*pi)/ 9; ‘angular frequency of Rotation 

AZ=abs(1.52e-8* (sin(wp*t+ps)+sin(wr*<t-+pe) jor 2; Sinusoidal Sun Area Variation 
A4=abs(1.52e-8* (sin(wp*t+ps)+sin(wr*t+pe)))/2; Sinusoidal Earth Area Variation 
A3=Al; Area of Probe Radiating 

A4=(D*L); $Area of Probe seen by Earth 

A5=A4; tArea of Probe seen by Earth reflected Sunlight 

Aé=Al; Area of Probe seen by Atmosphere 


Energy to/from the Probe: 


Solar=(alpha*AZ.*Ps); Energy from the Sun 
Planck=(epsilon*sigma*A3*T.%4); Planck radiation term for Probe 
EarthS=(alpha*A4*Psurf*gamma) ; Reflected Energy of Sun from Earth 
EarthP=(k*epsilon*A5* Pe); Long Wave Planck Radiation from the Earth 
AtmP=(epsilon*A6* Pa); $Long wave Planck Radiation from Atmosphere 


Tprobe=(Solar-Planck+EarthS+EarthP+AtmP+Al*h.*T) ./(h.*Al); Final Probe Temp 


$deltaT=.095; 
Convectieon—h. Al. “(lp rese— 1); Convective Heat Transfer from wire 


$semilogx (Tprobe-T,Altitude/1000); 

plot(h, Altitude/1000) 

$grid; 

$title('Probe Delta Temperature vs. Air Temperature'); 
ylabel('Altitude (km)'); 

xlabel('h'); 
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APPENDIX C - DATA ANALYSIS PROGRAM 


MATLAB Data Analysis program developed by Prof. Donald L. Walters: 


% 
$ 
% 
$ 


vizdata.m program to read the Sept 90 Anderson peak 


21 Apr 1997: "DLW 

set up constants 
fsample=1/1.25; 
£0=1/(2*pi); $fiicer 
We ue—2Z OCS 
AZ=NLLt/ 2: 
N=o2 


setup up file parameters 
fid=fopen('c:\Met data\Ap092716.txt','r'); 
cols=7; %$ # of columns in the data set 
col=5; % column # of temperature data 


read the file of arbitrary size 
(a, count) =fscanft(fid,*Sf\t"); 
Status=tclose(fid); 
rows=count/cols; 
a=reshape(a,cols, rows); 
a=transpose(a); 


extract the high speed temperature data column 
b=a( col); 
c=b(rows-(nfft-1): rows); tselect end of array 
c=detrend(c); 
G=1Ete (Cc); 
e-con] (a )°.=¢; 
ps=e(2:(n2+1)); tsubset of psd without DC term 


create frequency array 
£=26 ne 
f=f*fsample/nift; 
PlOGt =f, os); 


create RC high pass filter 
fat f7 £0). 2; 


ClGpem 7 (ltr) s $One pole fiiter 
fl Lr-—pilerec rer: Filter applied twice 
spectrai variance fitr * psd (with + and - frequencies) 


Ps2Z=cs- tier”; 
semilogy (it, ps2) 
Var=sqrt(2.*sum(ps2) ) 


time domain filter 
[ByAl=bureceril, 0. 16/7 .4 nigh”); $ create one pole RC 
OZ=filten(Sj,A, 5-52) )e % apply filter 
b3s=filter(8,A,b2); $ apply fiiter again 
var2=std(d3(100: rows) } 


Compute Cn2 
Z=a(:,2);3 % Altitude 
enzeb3..7bDs.* (7 9e=6" (ats 55). 7 (ate 4). 245.72; 


Lowpass smooth the data 
(De Gl=bDuUtEer( 2,205); 
Chie—< 1 perl Cnc) 
semi logxtenz, 2,0 6) 
axis({1.0e-20 1.0e-13 O 20000] ) 
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APPENDIX D- COEFFICIENT OF HEAT TRANSFER PROGRAM 


$Daniel J. Richardson 
$4 May 97 


%Coefficient of Heat Transfer Measurement 


%Define Constants 


Alpha=0.0042; Constant for Platinum Coated Tungsten Wire 
L=4.15e-3; Probe Wire Length 

D=4.7e-6; $Probe Wire Diameter 

R1=50; Resistance across Source 

RO=C.485; Probe Lead Resistance 

V1=2.0; Slnputc Voltage form Funetion Generator 


$Measured Voltages 


V2=0.2244; Voltage Measured at Low Point of Rise 
Ve=0223357; $Voltage Measured at Asymptotic High 
$Calculate Required Terms 

A=pi*D*L; tArea Of Probe Wire 

R2= (Rowan yyC aie oy $Probe lead Resistances 

R3= (Ras Vee V1—-V3S }eRO; + " 

dT=(R3-R2)/ (R2*Alpha) Differential Probe Temperature 


h=(V352Z) 7 (A-aT~ RS) 
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